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Introduction 

Ferrofluids  are  stable  colloidal  suspensions  o£  sub-domain  sized 
ferrite  particles  dispersed  in  a liquid  medium  by  a suitable  surfactant 
agent.  Ferrofluids  have  been  successfully  prepared  using  water,  hydro- 
carbons, esters,  diesters,  fluorocarbons,  and  even  liquid  mercury.  ; 

Since  their  recent  development,  considerable  Interest  in  the  possible 
application  of  ferrofluids  has  emerged.  A great  diversity  of  applica- 
tions has  been  proposed,  Including  propulsion  systems,  accelerometers, 
fluidic  elements,  "fluid  damping  components,  energy  converters,  and  mag- 
netic printers  .^Two  applications  showing  considerable  promise  are  ferro- 
fluid  rotary  shaft  seals  and  scrap  metal  separaters.  Rotary  shaft  seals 
have  been  commercially  available  for  several  years. 

Basic  research  on  the  magnetization  properties  and  the  bulkNresponse 
of  ferrofluids  to  magnetic  fields  has  been  conducted  to  meet  the  ne^^  of 
the  developing  technology  of  ferrofluids.  To  this  date,  research  on 
ferrofluids  has  been  fragmented  and  hampered  by  apparent  efforts  to  pro- 
tect proprietary  knowledge  concerning  the  constituents  of  these  fluids 
and  their  application.  Recently,  Dr.  B.  Berkovsky  has  organized  the 
first  meeting  devoted  specifically  to  ferrofluids  and  ferrofluid  appli- 
cations.* Many  things  about  ferrofluids  and  their  behavior  in  a magnetic 
field  are  not  yet  understood,  and  this  lack  of  knowledge  presently 
hampers  development  of  the  technology. 

^ From  the  fluid  mechanical  point  of  view,  ferrofluids  are  novel  because 
they  can  Interact  with  a magnetic  field  to  produce  a controllable  body  force 


International  Advanced  Course  and  Workshop  on  the  Thermomechanics 
of  Magnetic  Fluids,  sponsored  by  UNESCO,  Oct.  3-7,  1977,  Udine,  Italy. 
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on  the  fluid,  a body  force  significant  with  respect  to  terrestrial 
gravity.  From  the  basic  physical  point  of  view,  ferrofluids  are  inter- 
esting because  of  the  mechanisms  which  are  involved  in  the  transformation 
of  the  forces  on  individual  ferrite  particles  to  the  bulk  of  the  liquid 
carrier.  Historically,  two  areas  relevant  to  ferrofluids  have  received 
considerable  attention.  These  are  the  study  of  fluids  with  weak  magnetic 
properties,  and  the  study  of  solids  with  strong  magnetic  properties.  The 
wedding  of  these  two  fields  has  become  possible  and  necessary  through  the 
development  of  ferrofluids.  This  final  report  summarizes  three  years  of 
research  funded  by  the  Office  of  Naval  Research  under  contracts 
y/N00014-67-A-0299-0020  and  N0001A-76-C-0250.  .r>  The  researcli  program  was 
divided  into  studies  of  the  magnetization  proper ties  .^principally  investi- 
gated by  D.  A.  Krueger, -^and  the  electromechanics  and  applications^  princi- 
pally investigated  by  T.  B.  Jones,  both  of  Colorado  State  University.  An 
outline  of  the  essential  features  and  sigif leant  results  of  the  research 
is  provided  below  for  convenience  in  reviewing  the  program.  A more  de- 
tailed summary  follows  in  the  next  several  sections  of  this  report. 
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Magnetization  Studies 

The  principal  result  of  this  effort  has  been  the  discovery  of  a 
reversible,  magnetic  field  induced  agglomeration  effect  which  occurs  in 
water-base,  ester-base,  and  hydrocarbon-base  ferrofluids.  This  agglom- 
eration is  observed  in  constant  and  time-varying  magnetic  fields  as  low 
as  5 Oe.  It  occurs  whether  the  field  is  uniform  or  non-uniform.  Measure- 
ments Indicate  that  the  agglomerates  increase  to  sizes  involving  up  to  10® 

it 

Other  ferrofluids  were  not  examined  due  to  their  high  cost. 


Individual  ferrite  particles,  then  settle  out  due  to  gravitational  forces. 
The  agglomeratea  break  up  if  the  magnetic  field  la  rotiioved  and  gentle 
mixing  restores  the  ferrofluld  to  He  original  condition.  The  effect  le 
much  more  pronounced  in  water-base  ferrofluld  than  in  hydrocarbons.  Sub- 
stantial modeling  efforts  have  been  made,  and  some  of  this  effort  is 
reviewed  in  section  11. 


Electromechanics  and  Applications 


Simple  experiments  with  ferrofluids  to  determine  their  electro- 
mechanical response  to  magnetic  fields  have  been  conducted.  Principal 
among  these  is  the  Quincke  experiment,  in  which  the  height  of  rise  of  a 
column  of  ferrofluld  in  a magnetic  field  is  measured.  These  measurement.s 
are  found  to  be  consistent  with  magnetization  measurements  obtained  using 
a vibrating  sample  magnetometer.  Some  evidence  for  agglomeration  is 
noted  if  the  ferrofluld  Is  not  kept  mixed  throughout  the  measurements. 

A principal  conclusion  of  these  experiments  is  that  modeling  the  electro 
mechanics  of  ferrofluids  is  successful  only  if  the  magnetic  non-linearity 
Is  properly  accounted. 

A static  ferrofluld  plug  has  been  tested,  Mea.surements  made  with 
these  seals  are  consistent  with  a simple  ferrofluld  seal  model.  The 
mechanism  of  failure  for  a fluid/ fluid  seal  has  been  observed.  An  ex- 


tended series  of  experiments  with  a rotary  ferrofluld  seal  has  provided 
the  chance  to  gain  insight  concerning  the  important  seal  design  param- 
eters, including  ferrofluld  seal  inventory,  magnetic  flux,  gap  spacing, 
and  rotational  speed.  Thet-e  results  are  detailed  in  section  HI. 
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II.  Magnetization  Studies 

The  utility  of  and  intrinsic  interest  in  ferrofluids  arises  because 
these  fluids  have  magnetic  properties  as  well  as  fluid  properties.  We 
have  made  measurements  of  basic  magnetic  properties  for  several  fluids 
subjected  to  a variety  of  conditions,  and  we  have  used  magnetization 
measurements  as  a tool  to  investigate  agglomeration  phenomena.  Using 
our  newly  developed  experimental  technique,  we  have  observed  agglomera- 
tion in  concentrated  ferrofluids  like  those  used  in  devices.  Previous 
observations  of  agglomeration  were  made  in  very  dilute  ferrofluids  which 
are  poorly  suited  for  devices.  In  addition  our  observations  were  on 
fluids  with  measured  intrinsic  magnetic,  properties.  This  is  in  contrast 
to  earlier  experiments  and  aids  greatly  in  our  theoretical  understanding. 
A.  General  Propertle-s 

Using  a PAR  vibrating  sample  magnetometer,  we  measured  the  magneti- 
zation, M,  as  a function  of  applied  magnetic  field,  H,  for  various  combi- 
nations of  fluid  base  (water,  hydrocarbon,  and  ester),  storage  temperature 
(1“C,  25'’C,  and  SO'C),  and  magnetic  field  while  stored  (zero  field,  380 
Oe  uniform  field  and  a 500  Oe/cm  gradient  field).  After  about  6 months 
storage  the  magnetizations,  in  almost  all  cases,  had  dropped  by  about  5%. 
This  is  possibly  due  to  the  formation  of  agglomerates  which  have  some 
flux  closure  and  thus  do  not  respond  to  the  externally  applied  field. 

Osmotic  pressure  measurements  should  give  information  on  the  number 
density  of  the  magnetite  particles  in  the  colloid.  We  found  an  irrepro- 
ducibility  which  might  also  reflect  an  agglomeration  of  the  particles. 

An  electron  microscope  was  used  to  measure  the  size  of  the  magnetite 


particles.  We  found  that  the  water-base  fluid  had  significantly  larger 
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magnetite  particles  than  did  the  eater  or  hydrocarbon  base  fluids.  The 
size  distributions  for  these  two  fluids  were  indistinguishable.  Using 
the  Langevin  theory  to  calculate  the  magnetization  as  a function  of 
applied  magnetic  field  we  found  fairly  close  agreement  with  our  measure- 
ments  discussed  above.  However,  the  average  of  the  volume,  <v>,  as 
calculated  from  the  electron  micrograph  measurements  are  roughly  a 
factor  of  two  larger  than  those  deduced  from  the  M vs.  H measurements. 
This  is  not  surprising  because  a few  very  large  particles  observed  in 
the  electron  micrographs  contribute  very  heavily  to  the  average  volume. 

In  addition  the  extrapolations  of  the  M vs.  H are  uncertain.  These 
measurements  are  reported  at  length  in  earlier  annual  reports  of  this 
project  [1]  and  in  Peterson  [2].  These  measurements  indicate  that 
agglomeration  of  the  magnetite  particles  is  occurring.  We  next  briefly 
review  a long  series  of  experiments  which  unequivocally  demonstrate  this 
agglomeration. 

B,  Agglomeration  - Experiments 

The  basic  experiment  which  we  devised  to  study  agglomeration  is 
basically  to  measure  the  gravitational  settling  of  the  agglomerates  in  a 
vertical  tube  11  cm  in  length  and  0.6  cm  in  diameter.  A schematic  of  the 
experiment  is  shown  in  Figure  2.1.  The  shift  in  the  resonant  frequency 
of  the  Colpitts  oscillator  circuit,  Af,  upon  insertion  of  the  tube  into 
the  coil  was  found  to  be  proportional  to  the  saturation  magnetization  of 
the  fluid  in  the  column  within  0.37  cm  of  the  center  of  the  coil.  This 
in  turn  in  proportional  to  the  amount  of  magnetite  per  unit  volume.  A 
typical  result  for  Af  as  a functlok.  jf  position  along  the  tube  is  shown 
in  Figure  2.2.  It  is  clear  that  magnetite  particles  are  moving  from 


Fig.  2.1  Experimental  arrangement  used  for  lu  situ  magnetization 
measurements  using  oscillator  circit. 
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the  top  of  the  tube  to  the  bottom.  Analyzing  this  as  a function  of  time 
shows  that  the  particles  are  traveling  at  a constant  velocity  which  we 
identify  as  the  terminal  velocity.  This  terminal  velocity  depends  upon 
the  uniform  magnetic  field  which  is  applied  to  the  fluid  as  shown  bji  the 
crosses  in  Figure  2.3.  The  solid  line  is  the  result  of  our  theoretical 
model  which  we  discuss  in  detail  in  Section  C. 

From  our  data,  several  features  of  the  agglomeration  are  evident. 

The  agglomerates  are  made  up  of  on  the  order  of  10^  magnetite  particles. 
In  fact,  our  technique  is  insensitive  to  agglomerates  of  less  than  abqut 
10^  particles.  The  time  for  an  agglomerate  to  form  is  about  40  sec.  The 
time  for  an  agglomerate  to  break  up  upon  removal  of  the  magnetic  field 
is  about  10  sec.  The  agglomeration  is  reversible.  From  electron  micro- 
graphs of  Che  fluid  taken  from  the  top  of  the  tvibe  and  from  the  bottom 
of  the  tube  after  agglomeration  has  taken  place,  we  find  that  all  sizes 
of  magnetite  particles  in  the  fluid  participate  in  Che  agglomerates. 
Agglomeration  and  gravitational  settling  can  change  the  concentration  of 
magnetite  particles  by  factors  of  two  within  an  hour.  Water-base  fluid 
exhibits  much  more  agglomeration  than  does  an  ester  or  hydrocarbon-base 
fluid.  The  fraction  of  the  magnetite  which  participates  in  agglomeration 
Increases  as  If  there  is  a threshold  for  agglomeration,  it  is  less 

than  5 Oe.  Agglomeration  occurs  for  both  D.C.  and  60  hz  A.C.  magnetic 
fields  and  for  H hori rental  and  vertical.  These  experiments  are  dis- 
cussed more  extensively  ir  Peterson  [ 2 ] and  Peterson  and  Krueger  [ 3 ] . 
Agglomeration  also  occurs  in  iionuniform  magnetic  fields. 


of  the 


C.  Agglomeration  - Theory 

We  have  calculated  the  terminal  velocity  of  the  agglomerates  as  a 
function  of  applied  magnetic  field.  The  results  were  shown  in  Figure 
2.3.  The  terminal  velocity  can  be  calculated  easily  once  the  shape  and 
volume,  v,  of  the  agglomerate  is  known  because  the  relevant  Reynolds 
number  is  much  leas  than  unity  so  the  flow  is  viscous  dominated.  We 
assume  the  shape  is  a prolate  spheroid  (i.e.  cigar-shaped)  aligned 

along  the.  vertical  magnetic  field.  The  aspect  ratio,  k,  is  defined  as 

/ 

the  ratio  pf  uhe  semi-minor  axis  to  semi-major  axis  and  is  thus  less 
than  or  equal  to  one  (equal  for  a sphere) . To  determine  k and  v we  need 
two  equations.  Firstly  we  minimize  the  total  energy  of  an  agglomerate 
in  a magnetic  field.  The  energy  has  two  contributions.  The  magnetic 
dipole-dipole  energy  tends  to  align  the  particles  into  long  chains  (i.e. 
makes  k very  small).  On  the  other  hand,  the  surface  tension  energy  tends 
to  favor  a spherical  agglomerate  shape  (k = 1) . The  second  equation  is  a 
force  balance  on  the  particle  on  the  trailing  tip  of  the  agglomerate. 

The  viscous  drag,  F^,  as  the  agglomerate  falls  through  the  medium  tends 
to  pull  the  last  particle  from  the  agglomerate.  On  the  other  hand,  the 
magnetic  dipole-dipole  force,  Fj^,  attracts  the  particle  to  the  agglomerate. 

For  a marginally  stable  agglomerate  we  will  have  = F . 

M D 

In  principle  these  two  equations  suffice  to  fix  the  volume,  v,  and 
aspect  ratio,  k,  for  an  agglomerate.  Before  this  can  be  done,  there  are 
several  auxiliary  problems  to  be  solved.  The  calculation  of  F^^  is  non- 
trivial because  of  the  large  shielding  by  the  large  agglomerate.  If  the 
agglomerate  v;ere  a sphere  of  radius  R and  the  particle  has  a radius  R^, 
then  Goren  [ A ] showed  that  the  drag  force  F^^  is  proportional  to 
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For  our  agglomerates  this  factor  is  roughly  10  so  the  drag  is 
“ 6 

roughly  10  times  the  Stokes  drag  on  an  isolated  particle.  For  non- 
spherlcal  agglomerates  which  we  almost  certainly  have,  no  calculations 
of  Fjj  were  available.  Extending  the  multipole-expanslon-collocation 
technique  of  Gluckman,  et.al.  [5],  we  found  that  for  relative  volumes 
of  less  than  10**  a reasonable  analytical  approximation  to  F^^  is 


F = 19.4it  n n R 

U O 


where  n is  the  viscosity  of  the  ferrofluid  and  u is  the  velocity  of  the 
fluid  at  a distance  from  the  end  of  an  isolated  agglomerate.  VJe  have 
extrapolated  this  expression  to  larger  size  agglomerates  where  the 
collocation  technique  falls  due  to  numerical  problems.  The  present 
Justification  of  this  extrapolation  is  based  on  a reflection  technique 
as  discussed  in  Happel  and  Brenner  [ 6 ] > Further  discussion  of  this 
work  is  given  in  a paper  submitted  for  publication  by  Liao  and  Krueger  [ 7 ] 

The  drag  force  is  proportional  to  the  ferrofluid  viscosity. 

Rosensweig  [ 8 ] has  shown  that  the  viscosity  increases  by  factors  of  3 
or  4 in  strong  magnetic  fields.  He  also  showed  that  the  viscosity  is 
shear  dependent.  We  have  fit  his  data  with  empirical  curves  and  have 
used  these  in  our  evaluation  of  the  drag  force. 

To  evaluate  the  magnetic  force,  F,,,  and  the  magnetic  dipole-dipole 
energy,  we  must  determine  the  arrangement  of  the  magnetite  particles 
within  an  agglomerate.  For  any  applied  field  Luttinger  and  Tisza  [9  ] 
have  calculated  the  magnetic  dipole-dipole  energy  at  absolute  zero  of 
temperature  for  particles  arranged  in  three  lattice  structures,  simple 
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cubio,  face-centered  cubic  and  body-centered  cubic.  We  have  shown  that, 
for  all  applied  fields,  the  simple  cubic  structure  has  the  lowest  energy. 
Thus  we  assume  a simple  cubic  structure.  Sullivan  and  Deutch  [10]  have 
analyzed  a simple  cubic  structure  of  magnetic  dipoles  at  zero  applied 
field  but  for  nontero  temperature.  We  have  essentially  interpolated 
from  these  two  theories  to  obtain  results  for  nonzero  temperatures  and 
nonzero  fields. 

The  result  of  all  this  is  a rather  complicated  set  of  five  highly 
nonlinear  equations  which  must  be  solved  by  iteration.  The  surface 
tension  coefficient  is  fixed  by  fitting  the  theory  to  experiment  at  the 
lowest  fleJ.d  point.  Despite  the  complexity,  the  results  as  shown  in 
Figure  2.3  are  qualitatively  correct  and  quantitatively  quite  good. 

In  Figures  2.4  and  2.5  we  show  the  results  for  the  aspect  ratio, 
k,  of  the  agglomerates  and  their  volume  as  a function  of  applied  magnetic 
field.  The  particles  arc  very  elongated  (k  “ lO”^)  for  quite  small 
fields  (-20  Oe) . The  number  of  particles  in  an  agglomerate  is  very 
large  (-10^®)  for  similar  fields.  As  expected  the  agglomerates  become 
larger  and  more  elongated  as  the  magnetic  field  is  increased.  A more 
detailed  discussion  of  the  theory  is  given  in  Liao  [11]  and  Liao  and 


Krueger  [12]. 

D.  Areas  for  Further  Investigation 


Perhaps  the  most  important  area  for  further  work  is  to  investigate 
the  effect  of  shear  stresses  on  agglomeration.  Agglomeration  in  a 
static  fluid  is  important  because  many  devices  such  as  seals  or  magnetic 
ink  jet  printers  will  be  turned  off  for  long  periods  of  time.  In  addi- 
tion the  effects  of  agglomeration  In  a dynamic  setting  are  potentially 
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Fig.  2.4  Aspect  ratio,  k,  of  the  agglomerates  vs.  applied  magnetic  field 


Important.  Effort  should  go  Into  deciding  if  there  Is  a threshold 
shear  stress  for  breaking  up  agglomerates.  Additional  effort  should 
go  Into  understanding  the  shear  and  magnetic  field  dependence  of  the 
viscosity  of  ferroflulds. 


III.  Electromechanics  Studies 

The  essential  feature  of  the  response  of  ferrofluid  to  a magnetic 
field  is  electromechanical.  It  is  attracted  into  regions  of  more  in- 
tense magnetic  field.  The  force  is  exerted  on  the  individual  permanent 
magnetic  dipoles  of  the  ferrite  particles  by  an  imposed  non-uniform 
magnetic  field.  This  force  is 


(3.1) 


where  = 4tt  ' 10  H/ra  is  the  permeability  of  free  space,  in^  is  the 
magnetic  moment  of  a ferrite  particle  and  H is  the  externally  imposed 
magnetic  field.  The  net  force  per  unit  volume  may  be  obtained  by 
summing  the  forces  on  all  the  individual  particles: 


F = p M • 7H, 
o 


(3.2) 


where  M is  the  magnetization  vector.  Eq.  (3.2)  neglects  certain  inter- 
particle  (dipole-dipole)  interactions  which  have  no  effect  if  one  only 
wishes  to  calculate  the  total  force  on  a mass  of  ferrofluid.  Another 
more  common  force  density,  attributed  principally  to  Korteweg  and 
Helmholtz,  and  derivable  from  energy  considerations  [13] is: 


F = -%H  Vp 


(3.3) 
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vhere  p is  the  magnetic  permeability  of  a linear  ferrofiuid.  Most 
ferrofluids  are  not  linear,  requiring  a more  generalized  derivation. 

The  general  result,  valid  for  non-linear  ferrofluids,  is  [14] 

F - - E Va,  , (3.4) 

i *^1  ^ 


(3.5) 

and  u is  the  generalized  (non-linear)  permeability  as  defined  by 

1 " p(a^,  a^.  . . . , H^)H.  (3.6) 

The  are  fluid  properties  such  as  density,  ferrite  particle  concen- 

tration, etc.  Note  that  in  Eq.  (3.3)  and  (3.4)  rate-dependent  and 
hysteretic  magnetization  effects  are  neglected  and  the  continuum  approxi- 
mation for  the  colloidal  suspension  is  assumed. 

A.  Height  of  Rise  Experiment 

To  gain  confidence  in  the  behavior  of  ferrofiuid,  a classical  experi- 
mental technique  attributed  to  Quincke  [15]  was  repeated  using  a variety 
of  types  of  ferrofluids.  The  method,  used  to  measure  the  magnetic 
properties  of  weakly  susceptible  liquids,  involves  the  accurate  measure- 
ment of  the  height  to  which  a column  of  ferrofiuid  will  rise  in  a glass 
tube  when  subjected  to  a uniform  transverse  magnetic  field.  See  Fig. 

3.1.  The  height  of  rise  h is  a function  of  the  applied  uniform  magnetic 


magnetic 


Fig.  3.1  Illustration  of  experimental  apparatus  for  determining  the 
magnetization  properties  or  ferrofluiu  sampleB  using  the 
Quincke  method- 
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field  H 


P z*  H 

h-^/  M(H*)dH',  (3.7) 

P8  Jq 

where  p la  the  density  of  the  ferrofluid  and  g " 9.81  m/d^ec^  is  the 
acceleration  due  to  gravity.  The  density  of  air  is  neglected  here. 
Representative  experimental  results  for  two  ferrofluid  samples  are 
shown  in  Fig.  3.2.  As  a check  on  these  data,  a "theoretical"  height  of 
rise  was  calculated  using  M(H)  data  obtained  from  vibrating  sample 
magnetometer  data  (see  section  II)  in  Eq.  (3.7).  The  agreement  is  well 
within  the  uncertainty  of  the  measurements. 

The  above  measurements  were  carefully  made  using  only  well-mixed 
ferrofluid  samples.  Long  exposure  of  the  ferrofluid  to  the  magnetic 
field  affected  the  height  of  rise  data.  This  effect  depended  on  the 
length  of  the  exposure.  It  is  supposed  that  this  effect  is  the  same  as 
that  discussed  in  section  II. B above  and  attributed  to  particle 
agglomerat ion . 

B.  Parametric  Surface  Wave  Dynamics 

Parametrically  induced  surface  perturbations  may  be  important  in 
the  free  surface  pumping  of  ferrofluids,  and  so  parametric  surface  in- 
stability was  investigated.  A fairly  considerable  background  of  work 
on  the  interfacial  dynamics  of  ferrofluids  existed  prior  to  the  start 
of  this  research.  However,  parametric  surface  waves  had  not  been 
studied.  The  experiment  is  shown  in  Fig.  3.3.  The  parametric  surface 
instability  exhibits  a distinct  threshold  which  is  dependent  on  the 
frequency  and  the  magnitude  of  the  magnetic  field.  Above  this  threshold. 





perturbations  spontaneously  arise  which  grow  to  a certain  amplitude  and 
take  the  form  of  a standing  wave  pattern  with  an  easily  discerned  wave- 
length. This  wavelength  also  depends  on  magnetic  field  frequency  and 
magnitude. 

The  theory  of  interfacial  parametric  ferrohydrodynamics  is  anulogouf, 
to  that  of  electrohydrodynamlca  Il6],  except  for  the  non-linearity  of 
ferroflulds.  Therefore,  an  analytical  model  for  the  M(H)  relation  was 
required.  The  superparamagnetic  theory,  generally  considered  to  be 
accurate  for  ferroflulds,  was  employed  using  a single  Langevin  function; 


M(H)  “ M [coth(aH)  - l/(aH)] 
8 


(3.8) 


Eq.  (3.8)  Ignores  the  particle  size,  distribution,  substituting  an  average 
particle  size  which  best  fits  the  measured  M(H)  curve.  This  approach  is 
as  much  a curve-fitting  as  a modeling  exercise. 

Perturbations  of  the  surface  are  assumed  to  be  spatially  periodic. 
Using  Eq.  (3.8)  and  balancing  forces  at  the  interface  gives  the  surface 
equation  of  motion.  This  equation  fits  the  canonical  form  of  Hill's 


equation.  At  given  magnetic  field  frequency  and  magnitude,  a certain 
band  of  wavelengths  is  unstable  and  one  wavelength  value  has  the  faste.st 
growth  rate.  Wavelengths  predicted  by  this  "most  unstable  mode"  method 
are  compared  to  the  experimentally  observed  standing  wave  pattern  in 
Fig.  3.4.  Given  the  approximations  used  in  the  theoretical  model,  the 
agreement  is  reasonable,  suggesting  that  the  "most  unstable  mode"  hypoth- 
esis accurately  describes  the  development  of  the  non-linear,  amplitude- 
limited  standing  wave  pattern. 


Excitation  Frequ 

Figure  3.4  Plot  of  experimental  data  and  most  unstac 
parametric  ferroliydrodynaraics  experiment. 
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C.  Ferrofluld  Seals 

i Interest  In  shaft  seals  for  power  plants  atul  propulsion  systems 

suggested  the  Investigation  of  ferrofluld  seals.  Presently,  a line  of 

^raall  (“  3 cm  diameter)  ferrofluld  shaft  seals  Is  commercially  available 

j and,  further,  larger  ferrofluld  shaft  seals  for  cryogenic  and  vacuum 

} systems  are  being  developed.  Very  little  literature  on  the  design  and 

j operation  of  ferrofluld  rotary  shaft  seals  exists.  The  study  reviewed 

I below  was  initiated  to  gain  understanding  of  how  these  seals  work,  to 

i 

j identify  important  design  parameters,  and  to  develop  a model  which  ex- 

plained the  observed  behavior. 

Static  Plugs 

First  experiments  were  conducted  with  static  seals  in  order  to  gain 
understanding  of  the  essential  features  of  ferrofluld  sealing.  Refer  to 
Fig.  3.5.  The  static  plug  Is  established  by  concentrating  the  magnetic 
field  along  a short  section  of  a glass  tube.  Ferroflui"’  is  injected  into 
this  region.  Then  the  plug  is  loaded  from  above  by  an  immiscible  liquid. 
Using  this  very  simple  arrangement,  the  important  features  of  static 
ferrofluld  plugs  could  be  observed.  Because  of  the  liquid  loading,  plug 
failure  is  not  catastrophic;  no  appreciable  ferrofluld  is  lost.  A 
failure  channel  as  illustrated  in  Fig.  3.6  forms  and  the  upper  liquid 
passes  through  the  seal  along  the  tube  90°  off  the  axis  of  the  imposed 
magnetic  field,  where  the  field  is  the  weakest.  These  leaks  can  usually 
be  stopped  by  Increasing  the  magnetic  field  strength. 

The  model  developed  to  describe  the  behavior  of  the  ferrofluld  plug 
is  based  on  the  so-called  modified  Bernoulli  equation  [17]: 

p + pv^/2  ~ j M(R')dH*  = constant. 


(3.9) 


Fig.  3.5  Ferrofluid  Plug  Experiment 


'"^FAILURE 

CHANNEL 


Fig.  3.6  Details  of  failure  of  ferrofluid  plug  when  maximum  loading 
capacity  is  slightly  exceeded.  Small  channel  is  formed  90 
away  from  poles  because  of  small  variations  in  magnetic 
field  within  the  plug. 
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where  p is  the  pressure  and  v Is  the  fluid  velocity  (zero  in  the  hydro- 
static limit  considered  here).  Using  Eq.  (3.9)  for  the  idealized  plug 
shown  in  Fig . 3.7,  we  have 


'l'  '”2-  l-oX  «(»')<>»' 
"1 


(3.10) 


where  and  are  the  hydrostatic  pressures  on  either  side  of  the  plug, 
and  Hj  and  are  the  essentially  tangential  magnetic  field  strength  mag- 
nitudes at  the  two  interfaces.  This  model  clearly  shows  how  the  plug  is 
displaced  as  it  is  loaded.  The  maximum  loading  occurs  when  the  inter- 
face on  the  high  pressure  side  is  located  at  the  maximum  field  intensity 
with  the  interface  on  the  low  pressure  side  in  the  region  of  zero  magnetic 
field: 


•M(H')dH'. 


(3.11) 


Strictly  speaking,  this  limit  is  not  achieved,  but  in  a properly  designed 
plug  it  is  closely  approached.  Note  that  tiot  be  approached  if 

the  inventory  of  ferrofluld  is  too  small.  The  experimental  static  plug 
burst  pressure  is  compared  to  the  prediction  of  Eq.  (3.11)  in  Fig.  3.8. 

As  shown  in  Fig.  3.8,  the  model  described  above  successfully  pre- 
dicts the  burst  pressure  of  static  seals.  This  success  is  deceptive, 
because  the  theory  does  not  account  for  interfacial  stability.  Under 
certain  conditions,  Eq,  (3.10)  will  predict  successful  sealing,  yet 
experimentally  the  seal  fails.  This  apparent  discrepancy  is  Invariably 
due  to  interfacial  instability  which  permits  the  fluid  being  sealed  to 


-27- 


glass  TUBE 


MAGNETIC  POLE 


.FERROFLUID 

PLUG 


MAGNETIC  FIELD  INTENSITY 


POSITION 


OPTIMUM  INVENTORY 


^ EXCESS  INVENTORY 

INSUFFICIENT  INVENTORY 


''  \ ' 
\ ' ' 


POSITION 


\.  y 


M (H) 


H2 

/ M(H)dH 

Hi 


Fig.  3.7  a)  Static  ferrofluld  plug  between  opposite  poles. 

b)  Sjnnmetric  field  distribution  versus  position. 

c)  Hydrostatic  loading  for  varying  amounts  of 
ferrofluid  versus  seal  center  position. 

d)  Typical  M vs.  H relationship  for  ferrofluid. 
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rig.  3.8  Maximum  hydrostatic  Insding  of  ferrofluid  plug  versus  external 
(transverse)  flux  density  for  190  gauss  hydrocarbon-base  and 
250  gauss  water-base  ferrofluid.  Maximum  inventory  is  employed 
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pass  through  the  ferrofluld  layer.  This  behavior  is  more  common  in 
liquid-liquid  sealing  where  the  low  interfacial  tension  permits  the 
Rayleigh- Taylor  instability  to  destabilize  the  ferrof luid-liquid  inter- 
i face.  This  instability  can  be  suppressed  if  the  uniform  magnetic  field 

I 1 ^ ^region  ia  kept  very  small,  so  that  the  interfaces  are  stabilized  by  the 

non-uniform  magnetic  field.  This  matter  is  discussed  in  some  detail  by 

I ! 

{ Perry  [18]. 

i 1 

i 1 

Dynamic  Rotary  Seals 

j A simple  version  of  a rotary  shaft  seal  employing  ferrofluld  is 

1 

shown  in  Fig.  3.9.  This  seal  lacks  the  multiple  stages  commonly  em- 
ployed in  commercial  units,  but  it  is  useful  as  an  experimental  tool. 
i.  dynamic  seal  is  in  many  ways  similar  to  a static  plug;  however, ,the 
rotational  fluid  dynamics  complicates  the  device  behavior.  Burst 
pressure  measurements  for  the  device  shown  in  Fig.  3.9  were  made  under 

\ 

a variety  of  conditions.  These  included  variable  gap  spacings,  variable 
magnetic  flux,  and  variable  speed.  Some  representative  data  are  shown 
j in  Fig.  3.10.  The  different  curves  are  for  different  ferrofluld  inven- 

' tory  values.  Note  that  the  seal  performance  is  dependent  on  inventory 

as  predicted  by  the  static  seal  model  (see  previous  section) . Also  note 
‘ that  rotation  effects  are  prevalent  only  at  lower  inventory.  The  static 

burst  pressure  for  very  low  inventory  is  found  to  be  a small  fraction  of 
the  value  for  the  dynamic  (rotating)  case.  Both  the  rotation-dependent 
effects  seem  to  be  related  to  ferrofluld  redistribution  caused  by  centri- 
fugal forces,  however,  this  hypothesis  is  difficult  to  verify  due  to  the 
I difficulty  in  observing  the  ferrofluiu  location  within  the  seal  during 

rotation.  Likewise,  the  temperature  of  the  ferrofluld  is  difficult  to 
monitor  situ. 


! 
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Fig.  3. 


9 Cross-section  of  experimental  ferrcfluid  seal  and  test 
fixture.  Seal  is  loaded  from  below  by  evacuating  region 
above  seal. 
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Attempts  to  develop  a realistic  theoretical  model  for  the  seal  are 
frustrated  by  the  complex  nature  of  the  problem.  For  the  device  tested, 
magnetic  field  leakage  due  to  pole  block  saturation  made  the  true  loca- 
tion of  the  seal  uncertain  under  load.  Further,  the  influence  of  the 
rotation  on  the  fe.rrofluid  location  and  dynamics  was  dependent  on  the 
ferrofluid  location.  If  leakage  is  neglected,  the  model  is  not  capable^ 
of  explaining  many  of  the  observed  phenomena.  Thus,  the  modeling  exer- 
cise hinges  on  solution  of  the  non-linear  problem  of  ferrofluid  in  a gap 
with  strong  leakage  fields,  variable  Inventory,  and  variable  rotation 
speed.  As  yet,  only  Bailey,  et.al.  [19]  and  Perry  [18]  have  attempted  the 
task  of  modeling  a dynamic  seal.  Bailey,  et.al.  [19]  modeled  a conven- 
tional type  of  seal  for  a large  diameter  shaft,  using  an  energy  method. 
Perry  used  the  surface  force  balance  method  based  on  the  ferrofluid  plug 
model  to  analyze  a ferrofluid  seal  with  a transverse  magnetic  field  and 
a non-magnetic  shaft.  Perry's  model,  though  very  crude,  predicted  the 
effect  of  rotation  on  seal  burst  pressure.  No  one  has  yet  developed  a 
solution  for  a conventional  ferrofluid  seal  (with  magnetized  shaft) 


which  predicts  the  interfacia]  profile  of  the  ferrofluid. 
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IV.  Summary  and  Conclusions 

Over  the  three  year  period  of  the  research  contract,  the  basic 
magnetization  properties  of  ferrofluids  have  been  studied.  Reversible 
agglomeration  effects  not  previously  reported  have  been  observed  in 
vater-base^  eater-base  and  hydrocarbon-base  ferrofluids.  This  agglomera- 
tion is  dependent  on  .he  applied  magnetic  field,  and  it  occurs  in 
uniform  and  non-uniform  fields  as  low  as  5 Oe.  From  the  measured 

_ 3 

settling  velocities  of  approximately  10  cm/ sec,  the  apparent  agglomer- 
ate sizes  have  been  inferred  In  the  range  of  10^  particles.  The  agglomer- 
ates are  elongated  in  the  direction  of  the  applied  magnetic  field. 
Theoretical  predictions  of  the  maximum  agglomerate  size  are  based  on  a 

{ hydrodynamic  model  which  permits  drag  calculation  plus  an  energy  minimi- 

! ■ ■ 

I 

i zatlon  technique  which  permits  determination  of  the  particle  size  and 

shape.  The  work  is  substantiated  by  Independent  observation  of 
agglomerate  particles  made  by  Hayes [20].  The  agglomeration  of  water- 
base  ferrofluld  has  potential  implications  in  scrap-metal  separation 
technology  where  large  volumes  of  magnetic  liquid  are  employed  to 
separate  out  non-ferrous  metals  using  strong  magnetic  fields.  The 
requirement  of  cheap  ferrofluld,  the  large  volumes  involved,  and  the 

j strong  magnetic  field  imply  that  proper  racognition  of  agglomeration 

and  its  effect  on  magnetization  in  ferrofluids  is  essential.  The  useful 
static  shelf  life  of  certain  ferrofluld  seal  and/or  bearing  components 
is  also  a cause  for  concern  if  the  ferrofluld  is  subject  to  agglomeration. 
All  the  ferrofluids  tested  showed  a complete  recovery  from  agglomeration 
upon  gentle  mixing,  so,  once  recognized,  agglomeration  problems  are 
likely  to  be  easily  solved  in  any  application. 
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Over  the  same  period,  primarily  basic  studies  of  the  electromechanics 
of  ferroflulds  have  been  conducted  in  order  to  test  the  applicability  of 
present  continuum  models  for  the  force  on  a ferrofluid.  Miller  and 
Resler  [ 21]  have  questioned  the  generally  accepted  models.  Mo  experiments 
similar  to  those  of  Miller  and  Realer  were  attempted  in  this  study;  the 
expexlroents  conducted  did  not  indicate  any  breakdown  of  present  continuum 
models,  nor  did  they  rule  out  the  possibility  that  the  discrete  nature  of 
magnetic  liquid  colloids  exhibit  non-continuum  behavior.  Certain  experi- 
ments, notably  the  height  of  rise  measurements  did  uncover  a time-dependent 
effect  on  magnetization,  probably  attributable  to  agglomeration  which 
affected  the  ferrofluid  response  to  strong  fields.  This  effect  in  itself 
is  not  a breakdown  of  the  continuum  model,  but  rather  a feature  of  the 
magnetization.  However,  the  distinct  possibility  exists  that  the  large 
agglomerate  particles  might  lead  to  a breakdown  of  the  continuum  model. 

1 This  is  because  the  agglomerates  are  up  to  10^  times  larger  than  the 

Individual  sub-domain  ferrite  particles  themselves. 

The  parametric  surface  wave  dynamics  have  been  studied,  and  are 

j 

fairly  well  understood  at  this  point.  The  potential  for  application  of 
this  work  to  either  synchronous  interfacial  pumping  of  ferrofluid  or 

i 

magnetic  liquid  jet  printing  has  net  yet  been  pursued. 

A model  for  a static  ferrofluid  plug  has  been  developed  and  tested 
by  experiment.  The  model  is  sufficiently  simple  and  general  to  describe 
all  the  essential  features  of  ferrofluid  sealing.  Intensisre  experiments 
have  been  conducted  with  a single-stage  rotary  shaft  seal.  The  in- 

I fluences  on  the  burst  pressure  of  such  device  parameters  as  ferrofluid 

I 

I ' Inventory,  gap  spacing,  magnetic  flux  and  rotation  have  been  examined 


Iililtl 
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carefully.  Some  of  the  features  of  the  observed  behavior  have  been 
found  to  be  consistent  with  the  static  plug  model  and  some  have  not. 

The  complexity  of  the  rotational  seal  problem  has  thus  far  prevented 
the  successful  modeling  of  the  rotary  seal.  The  successful  application 
of  ferrofluid  seals  to  larger  shaft  sizes  and  higher  surface  velocities, 
as  proposed  for  shaft  seals  In  cryogenic  power  system  applications,  pro- 
pulsion systems,  etc.  depends  on  Improved  understanding  of  the  basic 
mechanisms  involved.  Considerably  more  work  on  the  modeling  of  dynamic 
ferrofluid  seals  Is  warranted. 


I 
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VI.  Other  Work 

A.  International  Advanced  Course  and  Workshop  on  the  Thermomechanics 
of  Magnetic  Fluids 

In  response  to  a request  from  Dr.  B.  Berkovsky,  one  of  the  Investi- 
gators (T.  B.  Jones)  has  prepared  a series  of  three  lectures  to  be 
presented  during  October,  1977,  at  the  "International  Advanced  Course 
and  Workshop  on  the  Thermomechanics  of  Magnetic  Fluids".  This  course 
is  sponsored  by  UNESCO  and  is  being  held  at  the  International  Centre  for 
Mechanical  Sciences  in  Udine,  Italy.  The  lectures  to  be  presented  con- 
cern the  theoretical  treatment  of  problems  in  ferrofluids,  with  emphasis 
on  application  to  seals.  The  lecture  series,  entitled:  "Theory  and 
Application  of  Ferrofluld  Seals",  is  divided  Into  three  sections  which 
are  1)  "Force  density  formulations  for  magnetic  fluids",  2)  Auxiliary 
fo*  a calculation  methods  and  solution  of  simple  ferrohydrostatic 
problems",  and  3)  "Static  ferrofluld  seals  and  plugs".  The  abstract  is 
reproduced  below: 

Theory  and  Application  of  Ferrofluld  Seals 

These  lectures  on  the  theory  and  application  of  ferrofluld 
seals  are  divided  into  three  topics  including  the  formulation 
of  forces  in  ferrofluids,  solution  of  various  ferrohydrostatic 
equilibria,  and  consideration  of  static  ferrofluld  seals  and 
plugs.  The  important  ferrofluld  force  formulations  include 
tigM  ■ Vff,  and  certain  force  densities  derived  from  energy  con- 
siderations. Each  of  the  force  densities  has  associated  with 
it  an  appropriate  stress  tensor.  In  the  case  of  incompressible 
magnetic  fluids,  ambiguities  in  the  force  density  arise  which 
can  be  shown  to  be  associated  with  the  arbitrariness  of  hydro- 
static pressure.  Despite  the  diversity  of  these  body  force 
formulations,  they  can  all  be  reconciled  and  shown  to  be  con- 
sistent if  integrated  over  the  entire  volume  of  the  ponderable 
body  in  question.  Application  of  these  various  force  formu- 
lations and  demonstration  of  their  mutual  consistency  is  best 
shown  by  working  various  simple  examples  of  ferrohydrostatic 
equilibL'la.  - CejlLaiii  fj£  Liiese  geoiueLLieS  may  be  sCilveu  for 
ferrofluids  with  arbitrary,  non-linear  magnetization  charac- 
teristics. Carefully  working  these  problems  provides  a solid 
basis  for  understanding  the  electromechanical  response  of 
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magnetlzable  liquids  In  non-uniform  magnetic  fields.  The 
single  most  important  application  of  ferrofluids  right  now 
is  in  vacuum  and  pressure  sealing.  Static  plugs  and  rota- 
tional seals  are  Illustrative  examples  of  f errohydrostatic 
and  f errohydrodynamic  equilibria  which  are  amenable  to  solution 
by  a number  of  different  approaches,  some  more  convenient  to 
use  than  others.  Models  developed  to  analyze  the  behavior  of 
static  ferrofluid  seals  lead  the  way  to  the  design  of  seals 
with  performance  characteristics  suited  to  various  applications. 


B.  Bibliography  for  Ferrofluid  Research 

£.  A.  Peterson  [23]  compiled  a bibliography  of  176  references  on 
ferrofluid  research,  which  Perry  updated  in  1976 [18],  adding  26  new 
references.  This  bibliography  has  been  updated  once  again,  by  the  addi- 
tion of  26  new  papers.  The  bulk  of  these  papers  have  appeared  within 
the  past  15  years,  indicating  widespread  Interest  and  activity  in  ferro- 
fluid research. 
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6.  Curtis,  R.  A.,  "Second  Order  Effects  in  Fluidmagnetic  Buoyancy". 
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